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EXECUTIVE  SUMMARY 


The  purpose  of  this  work  is  to  develop  new  approximate  theories  for 
scattering  from  randomly  rough  surfaces  which  are  less  restrictive  than 
conventional  approximations.  The  approach  used  was  to  apply  the  method  of 
smoothing  to  the  Fourier  transformed  Magnetic  Field  Integral  Equation  to  obtain 
a  new  integral  equation  for  the  fluctuating  part  of  the  scattered  field.  This 
new  integral  equation  contained  a  new  Born  term  which  was  the  sum  of  the 
transformed  Kirchhoff  approximation  and  a  new  term  proportional  to  the 
transform  of  the  fluctuating  part  of  the  surface  propagator.  This  new  Born  term 
reduced  to  known  asymptotic  results  in  both  the  high  and  low  frequency  limits. 
Further  study  is  required  to  fully  understand  how  this  new  Born  term  increases 
the  range  of  validity  of  our  present  approximations.  However,  it  appears  that 
it  seems  to  remove  the  small  slope  approximation  contained  in  the  Rice 
approximation.  If  so,  this  result  will  greatly  help  in  the  understanding  of 
scattering  from  small  scale  waves  on  the  ocean  surface. 


Abstract 


This  report  develops  an  integral  equation  for  the  fluctuating  zero  mean 
part  of  the  field  scattered  by  a  randomly  rough  perfectly  conducting  surface. 
The  development  starts  with  Magnetic  Field  Integral  Equation  for  the  current  and 
then  Fourier  transforms  this  equation  to  obtain  an  integral  equation  fo^  the 
scattered  field  in  transform  space.  The  method  of  smoothing  is  then  used  to 
develop  an  integral  equation  for  the  fluctuating  part  of  the  scattered  field. 
The  Born  term  in  this  integral  equation  of  the  second  kind  contains  much  more 
information  than  the  classical  Kirchhoff  approximation.  It  is  shown  to  reduce 
to  the  Kirchhoff  approximation  for  gently  sloping  surfaces  with  large  rms 
roughness  and  the  Rice  perturbation  result  for  small  surface  roughness.  This 
new  Born  term  contains  factors  which  appear  to  have  the  potential  capability  for 
going  beyond  existing  approximations  but  further  analytical  and  numerical  work 
is  required  to  quantify  this  hypothesis. 
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1.0  INTRODUCTION 


Scattering  of  electromagnetic  energy  by  randomly  rough  surfaces  such  as 

the  ocean  surface  has  many  important  applications.  All  of  these  axe  based  on  the 

link  between  the  scattered  field  behavior  and  the  character  of  the  rough 

surface.  In  its  most  basic  form,  this  link  is  established  through  the  physics 

represented  in  Figure  1.1.  The  surface  is  illuminated  by  an  incident 

electromagnetic  field  .  This  field  induces  a  surface  current  Jg  on  the 

-♦ 

surface  which  reradiates  the  scattered  field  E  .  The  total  field  above  the 

s 

*♦  ■>  -* 
surface  is  E  +  E .  .  The  knowns  in  this  problem  are  the  incident  field  E.  and  the 

statistics  of  the  surface.  It  is  important  to  note  that  it  is  the  statistical 

characteristics  of  the  surface  which  are  known  and  not  any  or  all  realizations 

of  the  surface.  The  purpose  of  reviewing  these  basics  is  to  point  out  that  the 

best  one  can  hope  to  solve  for  are  the  statistical  characteristics  of  the 

scattered  field. 

1. 1  Statement  of  the  Problem 

In  order  to  calculate  these  characteristics  of  the  scattered  field,  it  is 
first  necessary  to  treat  the  problem  as  deterministic.  This  is  usually  done  by 
setting  up  an  integral  equation  for  the  surface  current  induced  by  the  incident 
field  in  terms  of  the  incident  field  and  the  shape  of  the  surface.  Assuming  that 
th:s  integral  equation  can  be  solved,  the  resulting  current  weighted  by  a  phase 
factor  related  to  the  surface  height  £(x,y)  is  then  Fourier  transformed  with 
respect  to  the  x  and  y  surface  coordinates  to  yield  the  far  field  portion  of  the 
scattered  field.  This  part  of  the  scattered  field  is  adequate  when  the  field  is 
measured  sufficiently  far  from  the  surface;  a  condition  that  is  usually  always 
met  in  practice.  With  this  knowledge  of  the  scattered  field  and  its  dependence 
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SURFACE  Z  =  f(X,y) 


Figure  1.1.  Representation  of  an  incident  field  EF  impinging  on  a  surface 

defined  by  z=((x,y)  and  inducing  a  current  J  on  the  surface  which  reradiates 

s 

the  scattered  field  E  . 

s 
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on  the  surface  shape,  fundamental  theorems  of  probability  theory  can  be  used,  in 

principle,  to  compute  any  statistical  moment  of  the  scattered  field.  Of  course, 

what  is  really  desired  is  the  two— point  probability  density  function  (pdf)  of 

the  scattered  field  or  even  the  less  complicated  single  point  pdf.  The  former 

provides  estimates  of  how  far  one  must  move  in  space  before  the  scattered  field 

decorrelates.  Unfortunately,  neither  of  these  pdfs  are  very  easy  to  calculate 

because  of  the  complicated  manner  in  which  a  pdf  is  related  to  the  functional 

dependence  of  Egon  the  surface  characteristics.  Thus,  even  with  detailed 
-♦ 

knowledge  of  how  E  depends  on  the  shape  of  the  surface,  it  may  be  practical  to 
s 

calculate  only  the  first  few  statistical  moments  of  the  scattered  field 
including  its  temporal  and  spatial  decorrelation  properties. 

The  description  given  above  concerning  calculating  the  statistical 
properties  of  the  scattered  field  is  predicated  on  the  very  important  caveat 
that  the  integral  equation  for  the  induced  surface  current  can  be  solved.  This, 
of  course,  is  not  an  easy  task  and,  in  fact,  is  the  major  difficulty  in  this 
entire  problem.  If  the  surface  current  density  can  be  determined  then  all  other 
quantities  of  interest  can  be  obtained  via  straightforward  numerical 
integrations.  It  is  essential  however  that  the  integral  equation  for  the 
current  be  solved  to  a  known  level  of  accuracy. 

1.2  Possible  Solutions 

The  number  of  options  available  to  solving  the  integral  equation  for  the 
current  are  rather  limited.  One  can  employ  approximate  techniques  which  are 
usually  related  to  high  or  low  frequency  asymptotic  methods.  For  example.  Brown 
[1.1]  has  shown  how  the  current  induced  on  a  perfectly  conducting  rough  surface 
behaves  in  the  high  frequency  limit  while  Holliday  [1.2]  obtained  the  behavior 
of  the  current  in  the  low  frequency  limit  from  the  same  integral  equation. 
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Although  it  has  not  been  carried  through  to  completion,  it  is  not  difficult  to 
see  how  these  two  asymptotic  analyses  could  be  appropriately  combined  to 
obtained  a  first  order  estimate  of  the  current  induced  on  a  surface  having  two 
scales  of  roughness.  However,  since  such  a  combination  will  only  produce  a 
first  order  approximation,  it  is  doubtful  that  it  will  generate  results  which 
are  signif icantly  different  from  other  two  scale  theories  such  as  in  Brown 
[1.3] .  There  has  been  some  recent  work  appear  on  modeling  waves  as  wedges  and 
estimating  the  resulting  scattering  via  wedge  diffraction  [1.4]  .  However,  this 
particular  scattering  mechanism  appears  to  be  important  only  near  grazing 
incidence  and  it  depends  upon  the  statistics  of  wedge  like  waves  which  are  still 
largely  unknown.  In  general,  it  appears  that  approximate  solutions  for  the 
induced  surface  current  based  on  asymptotic  methods  have  yielded  about  as  much 
information  and  insight  as  they  are  capable  of. 

The  other  approach  to  solving  the  integral  equation  for  the  current  is  a 
pure  numerical  procedure  wherein  one  uses  the  methods  of  moments  to  convert  the 
integral  equation  to  a  matrix  equation  comprising  discrete  values  of  the  current 
along  the  surface.  Taking  the  Fourier  transform  of  the  current  yields  the  field 
scattered  by  this  particular  realization  of  the  surface.  Using  Monte  Carlo 
methods  to  generate  other  realizations  of  the  surface  from  a  known  density  and 
spectrum,  one  can  repeat  this  process  until  enough  scattered  field  values  are 
available  to  construct  the  desired  statistical  moments.  While  this  is  a  very 
powerful  method,  it  is  computationally  intensive  and  even  with  super  computers 
it  appears  to  be  practical  for  two-dimensional  (2D)  surfaces  only.  Two 
dimensionally  rough  surfaces  are  ones  where  the  surface  height  depends  on  one 
variable  only;  this  is  like  a  corduroy  surface  in  which  the  period  and  height  are 
random  functions  of  position. 

While  the  2D  surface  computations  do  produce  some  useful  results  and 
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insight,  it  is  not  clear  how  much  of  this  translates  to  the  more  general  3D 
roughness  case.  Clearly,  there  is  no  accounting  for  some  very  important 
polarization  effects  such  as  depolarization  or  the  conversion,  by  scattering, 
from  one  polarization  state  to  another.  This  particular  phenomenon  is  very 
important  because  it  forms  the  basis  for  remotely  sensing  certain  types  of 
surface  features.  There  are,  in  general,  no  guidelines  for  translating  2D 
scattering  results  into  3D  interpretations  and  so  it  is  necessary  to  consider 
solving  the  3D  problem  via  numerical  means.  Unfortunately,  this  dimensional 
extension  is  not  trivial.  This  is  because  one  must  overcome  the  difficulty 
associated  with  the  inversion  of  very  large  three  dimensional  matrices. 
Researchers  working  with  scattering  by  finite  size  bodies  have  found  that  a 
purely  numerical  approach  is  not  practical  when  the  size  of  the  body  exceeds  a 
few  electromagnetic  wavelengths  [1.5]  .  This  result  does  not  bode  well  for  the 
extended  rough  surface  problem  where  the  surface  is  necessarily  very  large 
compared  to  the  electromagnetic  wavelength. 

Given  this  very  practical  limitation,  what  options  are  available  for 
solving  the  general  3D  roughness  problem?  One  option  is  the  so-called  "Full 
Wave"  approach  of  Bahar  [1.6].  This  method  i*s  one  in  which  the  mathematical 
approach  is  not  at  all  obvious,  the  basic  physics  of  the  solution  are  obscure, 
and  the  first  order  approximate  results  are  claimed  to  contain  significantly 
more  than  similar  results  obtained  by  more  obvious  means.  In  short,  the  "Full 
Wave"  method  appears  to  raise  more  questions  (as  to  its  capabilities  and 
limitations)  than  it  answers.  For  this  reason,  the  method,  in  its  present  state 
of  development,  does  not  appear  to  be  an  attractive  alternative.  Certainly  this 
situation  could  change  if  the  method  were  more  carefully  explained. 

Another  possible  approach  might  entail  a  numerical  solution  of  Maxwell's 
partial  differential  equations  and  associated  boundary  conditions  by  some  form 
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of  mesh  grid  scheme.  This  technique  is  attracting  attention  in  wave  propagation 
problems  where  it  makes  maximum  use  of  the  computational  capability  of  so  called 
super  computers.  However,  it  is  not  completely  clear  how  this  approach  will 
perform  when  there  is  a  single  arbitrarily  shaped  interface.  It  is  not  at  all 
obvious  that  it  is  superior  to  the  integral  equation  formulation  and  so  this 
approach  requires  much  more  study. 

The  only  alternative  that  appears  to  be  practical  is  the  development  of  a 
new  integral  equation  which  has  a  much  more  robust  approximate  solution.  That 
is,  while  an  exact  solution  of  this  alternate  integral  equation  may  be  no  more 
feasible  than  with  the  original  integral  equation,  its  approximate  solution  may 
be  much  more  accurate  over  a  wider  range  of  parameters  than  that  for  the  original 
equation. 

1.3  Thrust  of  This  Report 

The  purpose  of  this  report  is  to  develop  a  new  integral  equation  such  as 
described  above.  This  new  integral  equation  will  be  obtained  by  a  sequence  of 
very  straightforward  steps.  The  starting  point  will  be  the  Magnetic  Field 
Integral  Equation  (MFIE)  for  the  surface  current  induced  on  a  perfectly 
conducting  surface  by  an  incident  field.  The  next  step  in  the  process  is  to 
Fourier  transform  the  terms  in  this  integral  equation  to  obtain  an  integral 
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equation  for  the  scattered  field;  a  reasonably  well  known  k— space  integral 
equation  formulation  for  the  scattered  field. 

The  next  step  entails  the  application  of  the  method  of  smoothing  to  this 
integral  equation.  This  amounts  to  splitting  the  scattered  field  into  an 
average  value  and  a  zero  mean  fluctuating  part.  This  apparently  trivial 
decomposition  has  a  number  of  very  positive  benefits.  First,  it  permits  the 
derivation  of  an  integral  equation  for  the  fluctuating  scattered  field  which  is 


exactly  what  is  needed.  That  is,  there  is  no  need  for  the  average  scattered 
field  because  it  propagates  away  from  the  surface  in  the  specular  direction 
only.  The  second  benefit  is  that  all  the  effects  of  the  mean  scattered  field  and 
the  mean  propagator  can  be  summed  and  this  leads  to  a  much  more  robust  or 
informative  Born  term  in  this  integral  equation  than  in  the  original  equation 
for  the  total  (mean  plus  fluctuating  parts)  scattered  field.  It  is  exactly  this 
Born  term  which  is  the  desired  result. 
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2 . 0  DEVELOPMENT  OF  AN  ALTERNATE  INTEGRAL  EQUATION 

For  this  initial  work,  it  is  desirable  to  keep  the  scattering  problem  as 
simple  as  possible  without  obscuring  the  basic  elements  which  are  attributable 
to  surface  roughness.  To  accomplish  this  goal,  the  case  of  a  perfectly 
conducting  rough  interface  will  be  considered.  Techniques  developed  to  deal 
with  this  problem  can  also  be  applied  to  the  imperfectly  conducting  interface 
[2.1]  but  the  complication  is  that  both  electric  and  magnetic  currents  must  be 
considered.  In  view  of  the  very  lossy  nature  of  sea  water,  it  may  be  possible  to 
use  approximate  boundary  conditions  to  simplify  the  lossy  interface  problem 
even  further  (with  some  additional  work)  . 


2 . 1  The  Basic  Integral  Equation  for  the  Current 

For  an  interface  which  is  perfectly  conducting,  there  are  two  distinctly 
different  kinds  of  integral  equations  for  the  surface  current  [2.2]  .  The  first 
is  derived  from  the  continuity  of  the  tangential  electric  field  across  the 
interface  and  it  is  called  the  Electric  Field  Integral  Equation  (EFIE)  .  It  is  an 
integral  equation  of  the  first  kind  which  does  not  lend  itself  to  an  obvious 
approximate  solution.  The  second  integral  equation  for  the  current  is  based  on 
the  discontinuity  of  the  tangential  magnetic  field  across  the  interface  and  so 
it  is  called  the  Magnetic  Field  Integral  Equation  (MFIE) .  It  is  an  integral 
equation  of  the  second  kind  whose  first  order  approximate  solution  is  outside  of 
the  integral  term.  Its  specific  form  is  as  follows: 


Vr)  =  Jg  +  2nW*J  Js(ro)*Vog(|r-ro|)dSc 


(2.1) 


where  J  is  the  unknown  surface  current  induced  on  the  surface  S  ,  and 
s  o 
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(2.2) 


g(|r-ro|)  =  exp(— jko|r-r|)/47r|r— rQ|  . 


In  the  above  n(r)  is  the  unit  normal  to  the  surface  at  the  point  r  =  xx  +  yy  + 

*  -♦  "j  -* 

(( x,y)z  on  the  surface  while  H  (r)  is  the  incident  magnetic  field  at  the  point  r. 

The  symbols  £  and  £  signify  the  x  and  y— slopes  of  the  surface,  k  is  the 
x  y  *  o 

wavenumber  for  free  space  (=2r/Ao),  and  g  is  the  free  space  Green's  function. 

■*  i 

Although  the  term  J  is  frequently  called  the  Born  term,  when  one  sets 


-f  -+  -j 

J  »  J1 
s  s 


(2.4) 


this  is  known  in  electromagnetics  as  the  Kirchhoff  approximation. 

The  MFIE  is  an  interesting  integral  equation  because  all  indications 

point  to  the  fact  that  the  Kirchhoff  approximation  is  not,  in  general,  as 

accurate  as  one  might  expect  [2.3].  That  is,  for  surfaces  containing  many 

scales  of  roughness,  the  Kirchhoff  approximation  does  a  good  job  of  predicting 

the  scattered  power  about  the  specular  direction  but  not  in  other  directions. 

■*  2 

This  is  due  largely  to  the  very  strong  coherent  effect  of  the  phase  of  H  in  the 
specular  direction  and  the  concomitant  lack  of  sensitivity  to  the  surface 
structure  comparable  to  A  .  The  reason  this  is  mentioned  here  is  that  it 
provides  a  first  clue  as  to  what  am  alternate  integral  equation  is  going  to  have 
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to  do  to  improve  on  (2.5).  That  is,  such  an  improvement  is  going  to  have  to 
deemphasize  the  strong  coherent  effects  of  Jg  in  the  specular  direction  and 
amplify  the  effects  of  the  small  scale  surface  structure  in  the  off— specular 
directions.  As  will  be  shown,  this  is  exactly  what  the  alternate  integral 
equation  does . 

2.2  The  Integral  Equation  For  The  Scattered  Field 

An  integral  equation  for  the  scattered  field  can  be  developed  directly 
from  (2.1)  in  the  so  called  far  field  approximation.  That  is,  if  the  point  of 
observation  of  the  scattered  field  is  sufficiently  far  from  the  surface,  it  is 
easy  to  convert  (2.1)  into  an  equation  for  the  scattered  field.  One  first 
multiplies  (2.1)  by  the  factor 


U  ♦  <vtc)2  exp  ()ks  o 


where 


and 

k  =  k  cos  9 
so  s 
z 


where  9  is  the  angle  of  the  scattering  direction  relative  to  the  normal  to  the 
s 

mean  surf  ace  (see  F .  g.  2 . 1) .  Eqn.  (2 . 1)  can  then  be  written  as  follows ; 


J(r)  =  J1  (r)  +  2N(r) 


N(r)xJ* 


J(r  ) 
o 


G(r-ro) 


dr_ 


(2.5) 
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Figure  2.1.  Incidence  and  scattering  geometry 
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where 


and 


J (r)  =  7  1  +  (V^Q2  exp  (jk  Q  J  (r)  (2.6) 

Z 

Jx(r)  =  2N(r)  x  Hx(r)  exp  (jk  Q  (2.7) 

S 

z 

N(r)  *  —V  (  +  z  (2.8) 


G(r-ro) 


exp 


jks  «-v] 

z  J 


(2.9) 


The  integration  in  (2.5)  is  now  over  the  z=0  plane  rather  than  the  surface;  this 
results  from  the  conversion 


dSo  -7 1  *  <w2  dK 

0 


-¥ 

where  dr.  =  dx  dy  . 
t  o  •'o 


The  second  step  in  the  conversion  process  is  to  note  that  the  current  must 
be  normal  to  the  surface,  i.e.  N  •  J  =  0,  or 


(2.10) 


Substituting  this  relationship  in  the  right  hand  side  of  (2.5)  leads  to  the 

following  matrix  equation  for  J  and  J  • 

x  y 
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J(r)  =  ji(r)  + 


*t  *C-C0.VtC0JtO  d^t 


(2.11) 


where 


11 

)  »“> 

'V 

K  * 

X 

i — 

>> 

J 

7„„ 

7, , 

L  y. 

'yx 

yy  J 

(2.12) 


and 


'  ac0  5g  _  dc  dg  dg  ' 

r  -| 

x^ 

X 

II 

to 

cfxQ  dx  +  dy  dy  d^ 

•  exp 

Jks  (C-C0) 

L  z  J 

(2.12a) 


'dC 

<?C01 

dg 

1 

7  *  —2 

'xy 

W 

3*  •  exP 

z  J 

(2.12b) 


'd(  dg  dc  dg  dg ' 

r  -j 

7  =2 

’yy 

•  exp 

jks  (C-C0) 

L  z  J 

_  ^C0 ' 

dg 

r  -] 

vf* 

X 

II 

A 

<?x  dx. 

0 

-$r  exp 

j'ks  (^_^o) 

L  z  J 

(2.12c) 


(2 . 12d) 


The  argument  of  the  propagator  K  shows  explicitly  which  variables  it  depends 

->  ■» 

upon.  The  notation  r.  —  r  implies  a  dependence  on  x— x  and  y— y  ,  i.e.  a 

XX  0  0 

o 

dependence  on  the  difference  in  the  transverse  (to  z)  variables.  If  (2.11)  can 

be  solved  for  J  and  J  then  J  can  be  found  from  (2 . 10) . 
x  y  z 

The  last  step  in  the  process  of  converting  (2.1)  into  an  integral  equation 

-♦  -♦ 

comprises  multiplying  (2.11)  by  exp(jk  -r  )  and  Fourier  transforming  the  result 
with  respect  to  x  and  y.  The  result  is  shown  below 
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/ru 

K(kt>  —  kt  )  -e(kt  )  dkt 
o  o 


(2.13) 


where 


(kt)  =J*  J(r)  exp  (jkt  -rt)drt 


(2.14a) 


ii(kt)  =j* Ji(r)  exp  (jkt  -rt)drt 


(2.14b) 


K^kt’  =  JJK(rt  rtQ’  C  Co*  Vo’  V5 


•  exp  j(kt  •  rt  - kt  *rt  )  drtdrt  (2.14c) 

oo  o 


The  term  in  (2.14c)  is  obtained  in  the  Appendix.  If  (2.13)  can  be  solved  for 


—  -*  ■» 

e(k  )  then  e  (k  )  can  be  found  from  transforming  (2.10),  i.e. 

X  Z  X 

e  (k  )  =  V  C®  e(k  ) 
z  t  t  t 


(2.15) 


where  the  tilde  over  V+C  denotes  its  Fourier  transform  and  ®  is  the  symbol  for 
convolution. 

A  solution  of  (2.13)  and  (2.15)  is  sufficient  to  determine  the  far  zone 

-»  -ft  ~  ~  ~ 

scattered  field  E  because  it  is  related  toe  =  e  x  +  e  y  +  ez  via 
s  x  y J  z 


Es  =  JVo*(R) V  Ve(kt  } 


(2.16) 
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or,  using  a  vector  identity  for  the  double  cross  product, 


Es  =  "j  VoS(R) 


e(k  )  -  [k  -e(k  )]k 
ts  s  ts 


(2.17) 


In  the  above  relationships,  t)q  is  the  characteristic  impedance  of  free  space,  R 

is  the  distance  from  the  origin  on  the  surface  to  the  point  of  observation,  g  is 

-• 

the  free  space  Green1  s  function,  and  k  is 

s 


k  =  kk  =  k 
s  os  o 


sin  6  cos  <f>  x  +  sin  6 


which  is  shown  in  Figure  2.1.  Note  that  in  order  to  determine  e  from  (2.13)  and 

■» 

2.15)  ,  one  must  find  it  for  all  values  of  k^ :  however,  only  its  value  at  kt  =  kt 

s 

is  required  to  determine  the  scattered  field. 


2.4  Using  Smoothing  For  The  Fluctuating  Scattered  Field 

Although  (2.13)  is  not  a  well  known  result,  it  fails  to  provide  the  kind  of 

—  jL  .  *♦  —  2. 

new  information  in  the  Born  term,  i.e.  e  (k^)  ,  that  we  seek.  This  is  because  e 

is  just  the  transform  of  the  Kirchhoff  approximation.  Something  additional  must 
be  done  to  convert  (2.13)  into  an  integral  equation  which  contains  a  more 
informative  Born  term.  To  accomplish  this  task,  the  method  of  smoothing  (MoS) 
will  be  used.  This  has  previously  been  used  to  determine  the  average  scattered 
field  [2.4],  but  it  will  be  applied  here  to  find  the  fluctuating  scattered 
field. 

To  simplify  matters,  (2.13)  is  written  in  a  more  compact  form,  i.e. 
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(2.18) 


(V 

e(k  )  =ei(kt)  +  LK(k  ,-kt  )  •  e(k  ) 

o 

where  L  is  the  integral  operator 

L  =  (27r)~2J*  (*)dkt 

The  key  element  in  the  Method  of  Smoothing  is  to  split  e  and  all  of  the  other 
quantities  in  (2.18)  into  the  sum  of  an  average  part  and  a  zero  mean  fluctuating 
part,  e.g. 


e  =  <e>  +  6  e 


(2.19) 


where  <e>  is  the  average  part  and  8  e  is  the  fluctuating  part  and  such  that 


<6  e>  =  0 


Decomposing  all  of  the  other  functions  in  (2. 18)  in  a  similar  manner  leads  to  the 
following  result : 


<e>  +  6  e  =  <ex>  +  6  e1  +  L<K><e>  +  L<K>£  e 


+  L6K<e>  +  L  8  K  8  e 


(2.20) 


Applying  the  fluctuating  operator,  1— <  >,  to  both  sides  of  (2.20)  yields 


8  e  -  8  e1  +  I.<K  >8e  +  L  £K<e>  +  (l-P)L  5K  6e 


(2.20a) 
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where  P  3  <  >  is  the  averaging  operator.  There  are  some  important 

ru 

simplifications  that  occur  in  (2.20a).  From  (2.14c),  the  average  of  K  can  be 
written  as  follows: 

<K<it  =  |  J  <iu;t  -;v  v- Vo» 


exp  j(k 


rt_kt 


drt  drt 


(2.21) 


When  K  is  averaged,  all  of  the  dependence  upon  (,  £  ,  Vt(,  aind  Vt(o  disappears. 

ss  +  -♦ 

What  is  left  in  <K>  is  a  dependence  on  r  —  r  only;  this  is  a  consequence  of  the 

o 

form  of  K  and  the  assumed  statistical  homogeneity  of  the  surface.  Thus,  (2.21) 
can  be  written  as  follows; 


<K(k  ,  —  k.  )> 


)> 

o 


exp 


dr. 


(2.22) 


With 


Art  =  rt  "  rt 


(2.21)  can  be  rewritten  as 
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<K(kt>  —  kt  )>  =  j  j <K(Art)>  exp(jkt  •  Art) 


•  exp  j(kt-kt  )  •  rt  dArtdrt 


o  o 


(2.23) 


which  becomes 


=  ”♦  -¥ 


<K(kt,-kt  )>  =  (2?r)^  <5(kt -kt  )  <K(kt)> 


(2.24) 


where 


(2ir)/«(k. 


Using  this  result  gives 


exp  j(k  -k  )  •  rt  dr 

O  0  0 


—  —  p 

L<K>  6  e  =  <K(kt)>  •  j  <5  (kt  —  kt  )  6  e  (kt  )dkt 

J  0  o 


L<K>  Se  =  <K(k  )>  •  5  e  (kt) 


(2.25) 


From  earlier  work  with  the  smoothing  method  [2.4]  ,  it  is  known  that  <e>  can 


be  written  as 


<e(k.  )>  =  (2x)  a(k  )  <5  (k  —  k  ) 
to  sz  to  ti 


(2.26) 
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2.5  Discussion  of  The  Resulting  Integral  Equation 

(2.28)  is  an  integral  equation  of  the  second  kind  in  transform  space  for 
the  fluctuating  part  of  the  scattered  field.  In  order  to  solve  it,  one  must  know 
the  average  scattered  field,  i.e.  the  a(k  )  in  (2.26).  However,  as  shown  in 

5 

z 

-* 

[2.4]  ,  a(k  )  is  itself  the  solution  of  an  integral  equation.  The  generation  of 
sz 

separate  integral  equations  for  both  the  average  and  the  fluctuating  parts  of 
the  scattered  field  is  the  real  power  of  the  method  of  smoothing.  The  reason 
this  is  important  is  that  the  separate  integral  equations  usually  lead  to  more 
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robust  approximations  than  first  making  an  approximation  in  (2.13)  and  then 
taking  the  average  and  fluctuating  parts  of  the  approximation.  With  the  latter 
approach,  it  is  seldom  clear  how  good  the  resulting  average  and  fluctuating 
approximations  are.  With  the  smoothing  approach  and  the  resulting  separate 
integral  equations,  it  is  much  easier  to  estimate  the  validity  of  a  given 
approximation. 

Thera  is  nothing  in  the  form  of  (2.28)  to  suggest  that  it  is,  any  general, 
an  easier  to  solve  than  (2. 13)  .  However,  it  does  contain  a  more  informative  Born 
term  than  (2 . 13)  . 

2.5.1  Signif  icance  of  The  Born  Term 

The  fluctuating  part  of  the  Born  term  in  (2.13)  is  the  fluctuating  part  of 
the  Kirchhoff  approximation.  As  discussed  earlier,  this  is  a  very  restrictive 
approximation.  The  Born  term  in  (2.28)  contains  more  information  than  just  the 
fluctuating  part  of  the  Kirchhoff  approximation.  The  factor  S  e  (k  )  with  k  = 

AJ 

k  is  the  Kirchhoff  terra.  The  factor  a(k  )  •  6  K(k  ,  —  k  )  with  k  =  k 

u  S  X  v  •  XX 

S  Z  1  s 

appears  to  bring  about  some  accounting  for  the  small  scale  structure  on  the 

surface.  In  fact,  when  it  is  combined  with  6  e  (k^  )  ,  the  sum  leads  to  the  Rice 

s 

2  2 

approximation  to  order  kQ  <(  >■  This  follows  from  the  form  of  the  sum,  the  fact 
that 


t\J  f\J 

$  K(k  ,  -k.  )  -  K (k  ,  -k,  ) 


for  small  surface  roughness,  and  the  work  of  Holliday  [2.5]  .  Thus,  the  Born  term 
in  (2.28)  contains  the  appropriate  behavior  for  small  surface  roughness. 
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As  the  roughness  becomes  large,  the  second  term  will  become  small  due  to 

-♦ 

the  decrease  of  the  mean  scattered  field,  i.e.  a(k  )  decays  almost 

z 

exponentially  with  increasing  roughness.  Furthermore,  if  the  slopes  are 

= 

sufficiently  small  then  <K(kt  )>  should  also  be  small.  Thus,  for  a  surface 

s 

2  2 

having  large  roughness  height,  <£  >,  but  small  slopes,  <(?  ()  >,  the  Born  term 
in  (2.28)  reproduces  the  proper  Kirchhoff  limit. 

For  surface  conditions  other  than  the  above,  it  is  necessary  to  study  the 

ru  aj 

behavior  of  both  <K>  and  8  K  because  they  will  obviously  bear  on  the  behavior  of 

AJ 

the  Born  term  in  (2.28).  The  factor  involving  <K>  appears  to  be  particularly 
important  because  it  involves  a  degree  of  apparent  multiple  scattering  on  the 
surface.  If  one  traces  this  term  back  to  its  source,  it  comes  from  the 
interaction  of  the  current  on  one  part  of  the  surface  with  the  current  at  another 

part  through  the  effects  of  the  average  propagator  <K>.  While  this  appears  to  be 
multiple  scattering  relative  to  the  current,  (2.28)  shows  that  it  has  a  single 
scattering  effect  on  the  scattered  field.  This  is  the  first  time  that  a  term  of 
this  form  has  been  identified  in  rough  surface  scattering. 
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Appendix 

The  purpose  of  this  Appendix  is  to  obtain  the  Fourier  transform  of  the 
integral  term  in  (2.11).  The  Fourier  transform  of  the  integral  term  can  be 
written  as 


I  = 


r  )  •J(ro)exp(jkt-  r  >dr 
o  o 


(A— 1) 


where  we  have  been  less  specific  in  showing  the  argument  of  K  than  in  (2.11). 
=  -♦ 

Since  only  K  depends  on  r^,  (A— 1)  can  be  written  as 


-  jVv 


Ktq) 


dr. 


(A— 2) 


Using  Parseval's  theorem  in  two  dimensions,  this  integral  can  be  written  as 
follows; 


rt  ^  ‘  J<r0>  drt 
o  o 


>  -J<*t 


CA— 3) 


>\l 


which  is  the  desired  result  because  J  =  e. 


I 
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3 . 0  CONCLUSIONS  AND  RECOMMENDED  FUTURE  WORK 

The  purpose  of  this  work  was  to  derive  a  new  integral  equation  for  the 
average  scattered  field  which  would  contain  a  more  robust  Born  term.  This  goal 
has  been  achieved  in  the  form  of  (2.28).  While  the  new  integral  equation  does 
not,  in  general,  appear  to  be  any  more  soluable  than  previous  integral 
equations,  it  does  contain  a  much  more  informative  Born  term.  This  term 
reproduces  both  the  Kirchhoff  and  the  Rice  perturbation  limits  and  appears  to 
have  the  capability  to  extend  beyond  these  asymptotic  limits.  In  fact,  this  new 
Born  term  has  the  potential  to  expand  our  existing  approximate  results  to  the 
point  where  further  numerical  work  on  the  full  3— D  integral  equation  may  not  be 
necessary. 

In  order  to  fully  exploit  this  new  result,  further  analytical  and 
numerical  work  is  required.  The  analytical  work  is  needed  in  order  to 
understand  the  physics  behind  each  of  the  factors  in  the  Born  term  in  (2.28); 
that  is,  their  source  and  meaning.  Such  insight  is  essential  to  estimating  just 
how  far  the  Born  term  can  be  pushed  before  it  becomes  an  inadequate 
approximation. 

Numerical  work  is  needed  to  determine  the  influence  of  the  average  and 

fluctuating  parts  of  the  propagator  K.  Since  their  influence  represents  the  new 
parts  of  the  Born  term  in  (2.28),  this  numerical  work  is  essential.  As  a  first 
cut  at  the  numerics,  the  2— D  form  of  the  Born  terra  in  (2.28)  should  be  compared  to 
some  of  the  numerical  solutions  of  the  MFIE  which  have  been  obtained  in  2— D. 
This  would  provide  some  guidelines  as  to  when  the  new  Born  term  is  accurate.  The 
evaluation  of  this  term  in  3— D  is  more  complicated  and  will  require  more  time. 
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